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Abstract 

We show that B-model topological strings on local Calabi-Yau threefolds are large N 
duals of matrix models, which in the planar limit naturally give rise to special geometry. 
These matrix models directly compute F-terms in an associated M = 1 supersymmetric 
gauge theory, obtained by deforming Af = 2 theories by a superpotential term that can 
be directly identified with the potential of the matrix model. Moreover by tuning some 
of the parameters of the geometry in a double scaling limit we recover (p, q) conformal 
minimal models coupled to 2d gravity, thereby relating non-critical string theories to type 
II superstrings on Calabi-Yau backgrounds. 
June, 2002 



1. Introduction 

Large iV limits of U(N) gauge theories have been a source of inspiration in physics, 
ever since 't Hooft introduced the idea In particular the large iV limit of gauge theories 
should be equivalent to some kind of closed string theory. The first contact this idea had 
with string theory was in the context of non-critical bosonic strings described by c < 1 
conformal field theories coupled to two-dimensional gravity. It was found that by taking a 
"double scaling limit" of N x N matrix models, where one send N to infinity while at the 
same time going to some critical point, one ends up with non-critical bosonic strings 
This relation between gauge systems and strings was not exactly in the sense that 't Hooft 
originally suggested for the large iV expansion, as it involved a double scaling limit. In 
particular, before taking this limit the matrix model would not have a string dual, whereas 
according to 't Hooft's general idea one would have expected it to have. 

In one context this was remedied by a different kind of matrix model introduced 
by Kontsevich Q, where without taking a double scaling limit one finds an equivalence 
between a matrix model and non-critical string theory. In particular the amplitudes of 
the topological string observables introduced in || are directly computed by these matrix 
integrals. This duality was in the same spirit as 't Hooft's general idea and can be seen as 
a low-dimensional example of a holographic correspondence. 

More recently large iV dualities have come back in various forms. In the context of 
M-theory a large iV matrix formulation was advanced ||. Since here an unconventional 
large iV limit is involved, this again was not quite in the same spirit as 't Hooft's idea, 
much as the double scaling limit of matrix models in the context of non-critical strings is 
not. However, the AdS/CFT correspondence |7j is in the same spirit as 't Hooft's original 
proposal in that one did not have to take a particular limit to obtain an equivalence. 
Another example of such a strict large N duality is the relation between Chern- Simons 
gauge theory and A-model topological strings || where one also does not have to take any 
particular limit for the equivalence to hold. 

The main aim of this paper is to develop a mirror version of this last duality || . We 
find matrix models that are dual to B-model topological strings on Calabi-Yau threefolds. 
This is again in the same spirit as 't Hooft, as one does not take a double scaling limit. 
We will show in particular that the special geometry of Calabi-Yau threefolds that solves 
the B-model at tree level emerges naturally from the dynamics of the eigenvalues of the 
matrix model. 
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However, even though it is not required, one can also consider a double scaling limit 
of this setup and obtain a specific class of Calabi-Yau manifolds that are dual to double- 
scaled matrix models. In this sense we are enlarging the original equivalence of double 
scaling limits of matrix models with string theory to an equivalence of all matrix models 
with some kind of closed string theory, without any need to take a double scaling limit. 
In particular our result shows that studying strings on non-compact Calabi-Yau spaces 
provides a unifying approach to all matrix model descriptions. 

Furthermore, it turns out that one can embed these large iV dualities in the context 



of type IIB superstrings || , a relation that was further explored in [pI)|p!l]Jl^ , |13| , P^ , p!5| . In 
the context of this embedding one obtains a dictionary in which the planar limit of the 
matrix models is seen to compute superpotentials for certain Af = 1 supersymmetric gauge 
theories, where the potential of the matrix model gets mapped to the superpotential for 
an adjoint scalar of an M = 1 theory. 

The organization of this paper is as follows: In section 2 we propose the large iV con- 
jecture for topological strings with matrix models after reviewing the various geometrical 
ingredients. In section 3 we pass this conjecture through some highly non-trivial checks, 
and in particular check it at the planar limit. In section 4 we discuss some generalizations 
of this conjecture and its connections with non-critical bosonic strings coupled to gravity 
and the double scaling limit. We also discuss the meaning of the double scaling limit from 
the viewpoint of type IIB superstrings. 



2. Large N Topological String Conjectures 

In this section we formulate a general class of large N conjectures involving B-model 
topological strings and certain two-dimensional topological gauge theories. The idea is 
to consider the mirror of the large N conjecture of |§, which relates large iV Chern- 
Simons theory on S 3 with A-model topological strings on the resolved conifold. This 
duality involves Calabi-Yau threefold transitions where a 3-cycle with branes wrapped 
over it shrinks and a 2-cycle grows without any branes wrapped over it. A worldsheet 
derivation of this duality has been recently presented in ||16|| . Moreover this has been 
generalized to a large class of Calabi-Yau threefolds |]T7[ which has been further studied in 

leading to development of powerful methods to compute all loop A-model topological 
string amplitudes [^9[^0[. 
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As suggested in || the mirror of A-model transitions dualities should also exist, namely 
we can consider 2-cycles with branes wrapped over them, undergoing transitions where they 
shrink and are replaced by 3-cycles without any branes. A large class of such Calabi-Yau 



transitions were considered in [TO|,|I^,piOJH| in the context of embedding such B-model 



dualities in type IIB superstrings. However, the duality between gauge theory and the 
topological B-model itself, which is implicit in these works, has not been studied. The aim 
of this section is to elaborate on these topological B-model/large N gauge theory dualities. 

2.1. Geometry of the generalized conifold transition 

Instead of being general we consider a special class of such transitions, studied in (TIJ , 



and discuss its topological lift. The situation considered in |10[ involved a string theory 
realization of AT = 2 supersymmetric U(N) gauge theory, deformed to M = 1 theory by 
addition of a tree-level superpotential TrVF($), which we take to be a general polynomial 
of degree n + 1 of the adjoint field $. The Calabi-Yau geometry relevant for this was 
studied in |Tl| and corresponds to considering the blowup of the local threefold given as a 



hypersurface 

uv + y 2 + W'(x) 2 = 0. (2.1) 

The blow up takes place at the critical points of W, i.e. at W'(x) = 0. Such transitions 
lead to a geometry that contains n blown up P^s which are all in the same homology 
class. So in the resolved geometry we can find n isolated rational curves. 

More precisely, the resolved singularity can be obtained by starting with the bundle 
O(0) © 0(—2) over P 1 . This is the normal bundle to rational curve in K3 x C and 
corresponds to M = 2 supersymmetry. Let us denote the sections of the normal bundle 
by </>o, (pi- These sections are respectively 0- forms and 1-forms on the P 1 . The field <po 
corresponds to the adjoined valued Higgs field $ in the M = 2 SYM theory. 

Now the inclusion of the superpotential W should give n isolated P^s at the critical 
values W'((j)Q) = 0. This is achieved by the following transition function. If z and z' are 
the coordinates on the northern and southern hemispheres of the P 1 , then the resolution 
is given by relating the patches (z, <fio,<fii) and (z 1 , </>q, 4>i) by 

z'<j>' 1 = z<j> 1 + W'(<fo), (2.2) 
z' = 1/z. 
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To relate it to the geometry given by ( |2.1|) we make the identifications 



x = </>o, 
u = 2<t>' 1 , 

v = 2<p 1 , (2.3) 



l! 

y = i(2u-W'((t> )), 

Here we have introduced another variable u in terms of which the geometry would have 
been given by 

uv -Auj 2 + 4uW'(x) = 0. (2.4) 

This makes a natural connection with objects that will be introduced in the next section. 

If we now distribute N D5 branes wrapped over these S 2, s and filling the spacetime, 
this corresponds to a choice of the vacuum in the corresponding Af = 1 supersymmetric 
gauge theory, where the distribution of the branes among the n critical points corresponds 
to distributing the iV eigenvalues of $ among the n classical values. At large N gaugino 
condensation takes place and this leads to a geometric transition in which the S 2, s are all 



blown down and replaced by "blown up" S 3 's. This results in the geometry |10 



uv + y 2 + W' 2 (x) + f(x) = 0, (2.5) 

where f(x) is a polynomial of degree n — 1, whose precise coefficients depend on the 
distribution of the N eigenvalues among the n critical points. 



In |15| it has been shown that, if one choses W(x) to be of degree N + 1 and if all the 
critical points are equally occupied, the M — 2 Seiberg-Witten geometry can be recovered 
by considering the limit W eW as e — > 0. 

Note in particular the identification of x with the eigenvalues of This is rather 
natural from the equation (|2.1|) . In particular if W = then we have an A\ geometry 
formed as a hypersurface in (u,v,y) space, where we have wrapped a D5 brane around 
the blowup sphere. The transverse position of the D5 brane along the x direction would 
correspond to changing the vev for $. Having a non-trivial fibration of the A\ geometry 
over the x-plane dictated by ( |2.1| ) gives rise to the superpotential VK($), as reviewed above. 



We will be interested in the B-model topological string on the deformed geometry 
given by equation fl2.5|) . The genus zero prepotential is determined by period integrals of 
the holomorphic (3,0) form 

dx A du A dv 



^uv + W'(x) 2 + f(x)' 

This non-compact Calabi Yau has n compact three-cycles Ai that correspond to the n 
"blown-up" S 3 and there are n canonically conjugated cycles Bi that are non-compact 
and have the topology of a three-ball. The usual special geometry relations determine the 
tree-level prepotential J-o(Si) by the periods 



'At 



dS 4 . 



(2.6) 

fi. 



Here one should be careful in regulating the periods over the non-compact B cycles |T(| . 



The degree of the deformation f(x), which corresponds to normalizable deformations, is 
such that these periods can be sensibly defined. In particular their variations with respect 
to the variation of the coefficients of f(x) are cutoff independent. 

As explained in [[Tof the A and B cycles can be represented as S 2 fibrations over paths 



in the complex x-plane. After integrating the 3-form over these 5" 2 fibers we are left with 
the integrals over these curves of the 1-form 

r\ = ydx. 

Here y is determined by the hyperelliptic curve 

y 2 + W'(x) 2 + f(x) = 0. (2.7) 

We will refer to this curve as the spectral curve. It has n branch cuts that are the 
projections of the cycles Aj onto the x-plane. The cycles Bi are represented by half-lines 
that start at the cuts and run to some cut-off point x = A far away from the branch cuts. 



2.2. Lifting to topological string dualities 

We now consider lifting these dualities to topological strings. As discussed in |§ the 
key point is the observation in |2!| that topological strings computes superpotential terms 
of the corresponding gauge theory arising in type II superstrings. In particular the leading 
planar diagram computes superpotential terms involving the gaugino bilinear field on the 
gauge theory side, and the Af = 2 prepotential on the dual gravity side (with some vev's 
for auxiliary fields, corresponding to turning on fluxes). More generally the topological 
gravity in the presence of Calabi-Yau 3-fold is the B-model theory studied in |2l| and 
called "Kodaira-Spencer theory of gravity" . Thus all we need to do is to specify the gauge 
theory dual, which should be the gauge theory on the topological branes wrapping the 
P^s. This theory is the reduction of the holomorphic Chern-Simons theory studied in [^3 
from complex dimension three to complex dimension one, which we now turn to. 

First suppose that we have no superpotential, i.e. W(x) = 0. In this case we simply 
get the Ai geometry times the x-plane. We wrap N branes around the P 1 . In this case 
the normal directions to the P 1 correspond to the cotangent bundle and the trivial bundle 
C associated to x. Let us call the Higgs fields in these two direction respectively &i(z) 
for the cotangent direction and <&o(z) for the x-direction. The topological theory on the 
B-brane we obtain in this case is given by the action 



S=- [ Tr($i£> A $ 



where $i(z) is a U(N) adjoint valued (1, 0) form on P 1 , <&o(z) is an adjoint valued scalar, 
A(z) is a U(N) holomorphic (0, 1) form connection on P 1 , and Da = d + [A, —]. Here g s 
denotes the string coupling constant. If we turn on the Higgs fields thereby deforming the 
P 1 to a non-holomorphic curve C, this action computes the integral 

S= — [ O 

9s Jy 

with Y a 3-chain connecting P 1 and C. 

If we turn on the superpotential W by deforming the geometry, the above action 
changes. This has been studied in |H] by studying the Beltrami differential that deforms 
the complex structure in the appropriate way, leading (including a minor generalization) 
to 

S w = — I Trf$!D A $ + W($ )u). (2.8) 
9s Jpi \ / 



where W(3*o) corresponds to the superpotential and w is a (1, 1) form which can be taken 
to correspond to have unit volume on P 1 . 

A consistency check for this action is to note that the equations of motion follow- 
ing from the above action agree with the fact that the classical solutions correspond to 
holomorphic curves. In particular integrating out A gives 

[$o,$i] = 0, 

so that $o and $i commute, i.e. we can assume they are simultaneously diagonal. Varia- 
tion with respect to the eigenvalues of $i leads to 

= 0, 

which together with compactness of P 1 implies that $o is a constant. Variation with 
respect to $o gives 

= w'($o)u 

which, together with the fact that the integral of <9$i over P 1 has to be zero for non- 
singular $i, leads to 

W'[$ ) = = $i. 

Thus the classical vacua indeed are localized at points W'($o) = 0, which describe the 
positions of the n P 1 ^.. 

In fact the action (|2.8D and the resulting quantum theory is rather trivial. In particular 
$i also appears linearly and can be integrated out exactly, leading to the constraint = 
0, which leads to the statement that $o is a constant N x N matrix 

$o(z) = $ = cnst. 

Thus the full action just reduces to its last potential term, which after integration of u 
over P 1 leads to the matrix action 

S w ($) = —TrW($). (2.9) 

9s 

Thus we see that the partition function of the gauge system is equivalent to a simple matrix 
model! 
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2. 3. Operator formalism 

We can give another derivation for the action (|2.8| ) starting directly from the patching 
functions ( [2.2| ) that determine the blown-up geometry. We split the P 1 in two hemispheres 
connected by a long cylinder. We denote the fields on these two patches as $o, $1 respec- 
tively $Q, In the case W = we are simply dealing with a gauged chiral CFT given 
by an adjoined valued ^-7 system of spin (1,0). The partition function computes the 
number of holomorphic blocks and this is one on the two-sphere. (Here we are ignoring 
for a moment the factor coming from the volume of U(N).) In an operator formalism this 
partition is simply given by pairing the left and right vacuum 

Z = (0|0) = 1. 

Now we want to implement the deformation induced by W. From the transition function 
fl2.2|) we see that the fields are related in the following way (here we write the fields in 
coordinates on the cylinder so that factors of z and z' are absorbed) 

Now there is an obvious operator U that implements this transformation on the Hilbert 
space. If we define 

U = exp — * Tr W($ (z))dz, (2.10) 
9s J 

(recall that the operator <&o(z) is an holomorphic field, so the contour does not matter as 
long as it encircles the poles) then one easily verifies that 

= U^tU' 1 . 

Here one uses the fact that the fields $0 an d $1 are canonically conjugated 

$o(z)$iW ~ 3s . 

z — w 

Therefore in an Hamiltonian formalism the partition function should be given by inserting 
the transformation U between the left and right vacua 

Z=(0\U\0). 

This is the familiar way to implement changes in complex structure in the operator for- 
malism of conformal field theory. 

In equation ( |2.10| ) we have written the deformation of the action in terms of the 
contour integration or Wilson line § TtW(<&q). Alternatively, this can be written as a 
surface integral 

U = exp— [ Tr^($ )w, 
9s Jpi 

where the volume form ui has been localized to a band along the equator of the P 1 . But, 
as noted above, one can take any 2-form, as long as it integrates to 1 over the sphere. 
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2.4- Precise Conjecture 

We are now ready to state our conjecture in precise terms. Let J-'wjiOs) denote the 
partition function of topological B-model for the Calabi-Yau manifold given by 

uv+y 2 + W'(x) 2 + f(x) = 0, 

where W is a fixed polynomial of degree n + 1 in x and / is a polynomial of degree n — 1 
in x. Let Si denote the integral of the holomorphic 3-form over the i-th S 3 coming from 
the i-th critical point of W. The periods Si will vary as we vary the n coefficients of /. 
Inverting this map, given the variables Si we can find the coefficients of the polynomial 
f(x) compatible with these periods. 

On the gauge theory side we now consider the matrix model given by the action 

S w ($) = — TrW(Q). 
9s 

We expand this matrix model near the classical vacuum given by partitioning 

N = Nx + ... + N n , 

and by putting iVj eigenvalues of $ in the i-th critical point of W . (Here we use both 
stable and unstable critical points. In fact, since we work in the holomorphic context, this 
difference does not really make sense.) Let ^ r vK,Af i (fl r s) denote the free energy of this matrix 
model expanded near this classical vacuum. Then the claim is 

Fwj{9s) = Fw,Ni{9s) 

with the understanding that Nig s = Si and, as discussed above, the periods Si fix the 
coefficients in the polynomial f(x). 

A few comments are in order: The matrix model integral is over a holomorphic matrix 
i.e. we have integrals of the form J <i$ and not like J <i$<i$. This is a general issue 
in topological B-branes0. The world-volume actions, such as the holomorphic Chern- 
Simons action in six dimensions, is a holomorphic function of the field variables. The 
non-perturbative path-integral should therefore by defined by picking some appropriate 
contour in the complex field configuration space. Indeed in our two-dimensional example 

1 We thank E. Witten for a useful discussion on this point. 
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we ended up with a chiral CFT and the path-integral definition of such a theory is rather 
subtle. One often ends up defining it as a holomorphic square root of a non-chiral theory. 
For a perturbative computation the situation is much easier. One simply performs the 
Feynman diagrams as if one was dealing with a real field. Similarly in this case we can 
effectively treat the matrix $ as if it is a real, that is Hermitean, matrix. 

Secondly we should note that the matrix integral, in the particular vacuum we end 
up, has a prefactor of 



This comes from the fact that for this vacuum U(N\) x ■ ■ • x U(N n ) denotes the unbroken 
gauge group and we have to mod out by the corresponding volume of the constant gauge 
transformations. This piece gives, as discussed in |16[] the partition function of c = 1 at 
self-dual radius. In particular the genus answer will involve 



In embedding in type IIB superstring this leads to the gaugino superpotential W e ff = 
^2 i NidJ r o/dSi + aSi. (Note that within the type IIB context the parameters Si and Ni 
are independent.) This is a first check on our conjecture. We are now ready to test the 
above large N conjecture in more detail. 

3. Matrix models and special geometry 

We will first show how our conjecture can be proven in the planar limit using standard 
manipulations in matrix model technology. A useful reference is for example ||. We will 
see how the special geometry of Calabi-Yau's emerges naturally. 

3.1. Matrix integrals in the planar limit 

We consider the one-matrix integral over N x N Hermitean matrices $ 



with W($) a polynomial of degree n + 1. By diagonalizing the matrix $ such integrals 
can be reduced to an integral over the eigenvalues Ai, . . . , Xn- In this way we obtain 



l/Vol(?7(iVi) x ••• x U(N n )). 



^ = E 2 S * logSi 
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where the Vandermonde determinants 



A(A) = n(A,-A J )=det(Ar 1 ) 

i<j 

appear from the Jacobian picked up by the diagonalization process. After exponentiating 
this contribution the effective action for the eigenvalues is given by 

s(\) = i £ W ^ - 2 E lQ g( A * - 

9s i i<j 

In this way we end up with a collection of N variables Ai, . . . , Ajv in a potential W(X). 
These eigenvalues interact through the second term that is a consequence of integrating 
out the off-diagonal components. This term gives the famous Coulomb repulsion between 
the eigenvalues. Because of the effective Pauli exclusion principle induced by the Vander- 
mondes it makes the eigenvalues behave as fermions. From this action we see that the 
equation of motion satisfied by the eigenvalues is 

9s ~f. *i — <\j 

We will now take the limit N — > oo of this system while keeping fixed the 't Hooft 
coupling 

fi = g s N. 

In this standard large N limit we will have a continuum of eigenvalues and their density 

i 

becomes a continuous function on the real axis normalized to f p(\)d\ = 1. (In the 
following A will always denote a real variable, in contrast with the variable x that can be 
complex). The eigenvalues will fill a domain on the real axis. This domain might consist of 
several disconnected components known as cuts. In the case of more than one components 
one speaks of a multi-cut solution. In the present case they are at most n of these cuts. 
We will denote the corresponding intervals in the complex x plane 
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3.2. The spectral curve 

To further analyze the model it will be convenient (and standard practice) to introduce 
the trace of the resolvent of the matrix <E> 

i 

This resolvent plays an crucial role in matrix model technology. It also has an interesting 
physical interpretation. For example, it can be thought of as a loop operator. 

By multiplying the equation of motion ( |3.2[ ) by the factor l/(Xi — x) and summing 
over i one obtains the important relation (loop equation) 



u 2 {x) - ^u'(x) + -u(x)W'(x) + -^f(x) = (3.3) 



where the polynomial f{x) is of degree n — 1 and is given by 

/(x)= £ E mz£w. ,3.4, 

i 

In some sense the function f(x) determines through fl3.3|) the whole solution of the matrix 
integral. Since it is polynomial of degree n — 1 we only have to determine the n unknown 
coefficients. 

In the large iV limit the second term in ( |3.3|) can be ignored and the differential 
equation for u(x) becomes an algebraic equation 

u 2 (x) + -u(x)W'(x) + -^f(x) =0 (3.5) 

From this we see that the resolvent u(x) in general has a piece that can have branch cuts. 
This singular part is captured by the function y(x) that we define here as 



y(x) = 2/ju(x) + W'(x) 
1 

Xi - x 



=2 9s j:^—+w'(x). (3 - 6) 



In terms of the variables (x, y) the relation (|3.5|) associated to the matrix model now takes 
the form 

y 2 - W'(x) 2 + f(x) = 0. (3.7) 
12 



This has an interpretation as a hyperelliptic curve in the (x, y)-plane. We immediately 
recognize this curve as the curve ( |2.7[ ) that determined the periods and thereby the tree- 
level topological string on the deformed Calabi-Yau geometry. In fact, it is even more 



natural to identify directly the parametrization of the curve in terms of (x, u) in (|3.5|) with 
equation ( |2.4| )of section 2. From this we see that we should identify uj with the field z<&\. 

Note that the function y(x) naturally appears from the dynamics of the eigenvalues. 
It is given by the variation of the action (|3.1|) with respect to a particular eigenvalue A 



y( X )=9s W (3.8) 

In fact, from this we see that the quantity y(x) has an elegant physical interpretation. 
Being the derivative of the potential energy, it equals the force on an eigenvalue if it tries 
to go away from its stationary position and moves into the complex a;-plane. 

3. 3. Filling fractions and periods 

We now want to evaluate the matrix integral around a particular stationary point 
where particular fractions of the eigenvalues cluster around the different critical points. 
Around such a multi-cut configuration it makes sense to make a perturbative expansion 
of the matrix integral using large iV techniques. It is the contribution from one of these 
saddle points that we are after. 

Consider such a multi-cut solution. The filling fractions Ni/N, i.e. the relative number 
of eigenvalues around each critical point, are given by the integrals 

Ni/N= [ p(X)d\ 

J A, 

Here the A\ are the cuts in the complex a;-plane. 

Now an important standard result is that the density of eigenvalues is given by the 
jump of the resolvent u(x) across the cut 

p(A) = — (u(X + iO) - u(X - iO)) . 

Since only the singular piece contributes, we can also write this relation as 

P (X) = -^(y(X + iQ)-y(X-iO)). (3.9) 
13 



Therefore we can compute the fraction of eigenvalues in a specific cut by doing a contour 
integral around the cut. In this way we find that the quantities Si = g s Ni are given by 
the period integrals around the cut, that is the periods on the Riemann surface 

Si = ^l y(x)dx. (3.10) 

Here we make contact with the period integrals ( |2.6| ) as obtained in the topological B- 
model computation. This is the first half of the derivation of the genus zero part of our 
conjecture. 

In order to complete the derivation, we now need to compute the change in the free 
energy J-o(Si) if we vary the filling fractions Si by adding an eigenvalue to the cuts 

ASi = —AN,. 

9s 

This change in the action is given by the work done by the force F(x) acting on an 
eigenvalue if we move this eigenvalue from one branch to infinity. We have seen in ( |3.8p 
that this force is given by 

F(x) = —y(x). 
9s 

So the variation of the free energy is computed in terms of the action J F(x)dx, that is by 
integrating the one-form y(x)dx along one of the S-cycles to the cut-off point x = A. We 
therefore immediately find the special geometry relation 

dF 



dS,. 



y(x)dx (3-11) 



B, 



expressing the S-periods in terms of the A-periods through the free energy !Fo(S). Together 
( p.l(J| ) and ( |3.11[ ) give the precise match with the Calabi-Yau geometry. This concludes 
the derivation of the planar version of our conjecture relating the matrix model to the 
topological B-model on the deformed Calabi-Yau. 

As we have mentioned, for a given potential W(x) the function f(x) that deforms the 
singular Calabi-Yau and determines the solution of the matrix model can be expressed in 



terms of the periods Si and vice versa. In fact, using definition ( |3.4| ) we can give a useful 
relation valid in terms of the prepotential JF . If one parametrizes the potential as 



n+l 



W(x) = ^u k x k , 



A;=0 
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and the deformation as 

71-1 

f(x) = ^4^6 fc x fc , 

fc=0 

then by plugging definition ( |3.4| ) in the matrix integral we obtain 

n+l 

b k = (k + 2)u k+2 + ^2 3 u o 



dUj-k-2 



j=k+2 

Note that if we introduce Virasoro operators L k = J2j J u jd/duj + k this equation can be 
written as p3 | 

bk = (k + 2)u k +2 + k = 0, • • • , n - 1. 

3.4- Higher genus 

Can one extend the derivation of our conjecture to higher genera? One possibility 
might be to use the loop equations that are derived by taking the expectation value of 
expression ( |3.3[ ) for the resolvent uj(x) 

(u 2 (x)) + -(u(x)} W\x) + -^f(x) = 0, 



where f(x) is now given by the expectation value of expression fl3.4|) . These loop equations 
have been studied for multi-cut solutions, for example in j55| , and one can try to solve them 
order by order in the string coupling constant g s . In principle these equations give recursion 
relations that relate the higher genus amplitudes in terms of the tree-level free energy. It 
would be very interesting to see if these equations are directly related to the equations of 
motion for Kodaira-Spencer string field theory [p2 j. This is not completely unlikely since 



collective field theories for the eigenvalue densities are known to have a similar form f26 
and morally the connection between topological strings and matrix models should go along 
these lines. 

It might be interesting to also briefly discuss the case of the pure conifold, here given 
by the quadratic superpotential W(x) = x 2 . The large iV dual is from our point of view 
just the Gaussian matrix model 

1 



<i$ • e" 



Vol(?7(iV)) 

Since the integral is trivial, the only contribution comes from the normalization factor 

N-l 

n * 



voK.n.Y)) / ,_ ] 
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which has been shown to reproduce the all genus answer for the B-model on the conifold 
in the 1/N expansion [|I6| |. In this case the spectral curve is given by 



y 2 — x 2 + [i = 

and the eigenvalue density 

p(X) = v/A^ 

is Wigner's famous semi-circle distribution. In the eigenvalue basis the all genus answer is 



alternatively obtained by the method of orthogonal polynomials which indeed gives |27 

N-l 



fjldX,- A(A) 2 ■ e -il A ? ^ Y[ fc 1 

k=i 



3.5. Domain walls and eigenvalue tunneling 

There is an interesting relation that directly connects the D-branes in the type II 
string theory and the behaviour of the eigenvalues in the matrix models. In the "old days" 



it was pointed out by Shenker |28[ that the characteristic non-perturbative effect observed 
in matrix models was the tunneling of eigenvalues and this was an effect of strength l/g s . 
This remark in some sense anticipated the importance of D-branes. Here we can connect 
the two effects. 

In the type IIB theory on the resolved geometry we can consider D5-branes wrapped 
around an S 3 that interpolates between two S 2 . Such an object manifest itself as a domain 
wall in the four uncompactified spacetime dimensions. After the geometric transition such 
a D5-brane will connect two three-cycles. It will describe a process where one unit of RR 
flux is transported. That is, the RR flux in one S 3 is decreased by one unit and the flux 
in another S 3 is increased by one. Since the space-time superpotential is given by 

w *ff=Y^( N *w ( +aSi )' ( 3 - 12 ) 

i 

the tension of a domain wall transferring flux from the i-th to the j-th cycle is given by 

<9JF dF , , . , 
y\x)dx. 



dSi dS j J Bi 



We now recognize this as the instanton action in the matrix model of an eigenvalue tun- 
neling from the cut Aj to the cut Aj along the path Bij. 
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4. Double scaling limits and other further generalizations 

Given that we have found a natural stringy interpretation of ordinary matrix model, 
one could ask what is the meaning of the double scaling limit in the context of the old 
matrix model 0. Following that limit on the gravity side for the single matrix model leads 
to the Calabi-Yau geometry 

H = uv + y 2 + x 2m+1 + deformations = 

as corresponding to the (2, 2m + 1) bosonic minimal model coupled to two-dimensional 
gravity. The deformations correspond to the m observables of the (2, 2m + 1) model. In 
fact for generic deformations of this geometry there are m A-cycles and m B-cycles and 
we thus can choose m independent parameters to parametrize these deformations. The 
infinitesimal deformations which map this geometry to the deformations of the (2, 2m + 1) 
models are of the form 

771 

uv + y 2 + x]J(x - e,) 2 = 0, (4.1) 
i=i 

where e$ are related to the deformation of the (2, 2m + 1) model with primary fields. For 
example, the case of pure 2d gravity, i.e. the (2,3) model has only one observable and in 
that case we have e = p}l 2 . In the usual matrix model this is obtained by considering a 
quartic superpotential W(x). Let us take that to be an even function of x. Then there 
are three critical points, including one at x = 0. Upon deformation by f(x) the critical 
points can split. If we put all the iV eigenvalues at the well corresponding to x = then 
only the x = double point splits and the other two do not split. The double scaling 
limit corresponds to taking the limit where one of the double points reaches one pair of 
doubles and the other reaches the other pair. Taking the limit where the geometry is 



localized near one set of triple zeros we obtain the (2,3) local geometry given in ( f4.1|) . 
Embedding this kind of theory in type IIB strings gives exactly the kind of N — ► oo 
dualities proposed recently in MM which relate certain limits of gauge systems with type 



IIB strings on Calabi-Yau geometries without fluxes. 

Note that this map to topological B-model is in perfect accord with the fact that the 
bosonic strings have a hidden Af = 2 superconformal symmetry [|29[], which in this case 



we can identify with the superconformal theory on the above Calabi-Yau threefold. As a 
check note that, if we only turn on the cosmological constant, then the genus g answer for 
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topological B-model scales as the holomorphic threeform O to the power of 2 — 2(/ [22 . 
For the (2, 3) model since we have 

dudvdydx 5 /4 
° = dH ~^ ■ 

So we learn that the free energy of the B-model is given by an expansion of the form 

^ = E c ^ 23 "V f(2 " 29) 

g>0 

for some coefficients c g , in perfect accord with the expected answer. 

As a further generalization one can consider the topological lift of the models con- 
sidered in \TA\. These involve many U(N) gauge fields with bifundamental matter. In 



particular for the case of a linear chain oip — 1 gauge fields one obtains the dual geometry 
(before blow up) 

v 

+ H(y-W[(x)) = 0. 



p 



where the Wi(x) are related to the superpotentials of the corresponding gauge fields. Again 
a double scaling limit exists which lead to geometries of the form 

uv + y p + x q + deformations = 

This would correspond to the (p, q) bosonic minimal model coupled to gravity. There are 
^{p — l)(q — 1) A-cycles for deformations in the above equation, which correspond to the 
observables of the (p, q) minimal model. As discussed in |13|] in this case the reduction of 
holomorphic 3-form leads to p — 1 1-forms rji which naturally get identified with the p — 1 
eigenvalue densities in a (p— 1) matrix model. It would be interesting to study these model 
in more detail and in particular the corresponding multi-matrix model duals. 

4-1. Final remarks 

There are many indications that topological strings on a general target space might 
be described by some kind of integrable systems. This was originally shown for the c < 1 
topological string theories PO] . More recently in the context of A-model topological strings 



on non-compact Calabi-Yau this integrability was shown by dualizing to certain observables 
of Chern-Simons gauge theory on three- manifolds [pTQ|j20| ] . For A-models on compact target 



spaces evidence has been accumulating in the mathematical literature on Gromov-Witten 
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invariants (see for example [31, 3^]). If we manage to formulate a matrix model dual of 



the topological string this integrability is in some sense manifest. Both the finite iV and 
the large iV matrix models are well-known to give tau-functions of the KP and Toda 
hierarchies [[3|,|27|- This integrability of matrix models was the underlying reason that non- 
critical strings with c < 1 were exactly solvable. Our results indicate that for a large class 
of non-compact Calabi-Yau manifolds this integrability is present. 

Furthermore one should also like to be able to include gravitational descendents within 
the topological string model. In terms of the B-model we expect that these are given by 
non-normalizable deformations of the complex structure. 

A related issue is the description of the c = 1 string, which is equivalent to topological 
strings on the conifold geometry |33|]. In a forthcoming paper |34[] we will describe how 



tachyon scattering processes are indeed reproduced in the conifold string theory and can 
be described by a large N dual gauge system, making contact with |35j and the recent 
work 
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